Dispersive microwave bifurcation of a superconducting resonator cavity incorporating 

a Josephson junction 
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We have observed the dynamical bistability of a microwave superconducting Fabry-Perot cavity 
incorporating a non-linear element in the form of Josephson tunnel junction. The effect, which is 
the analog of optical bistability, manifests itself in the transmission and reflection characteristics of 
the cavity and is governed by a competition between the wave amplitude dependence of the resonant 
frequency and the finite residence time of the field energy inside the cavity. This finite residence 
time is solely due to extraction of radiation from the cavity by the measurement process. The tight 
quantitative agreement with a simple model based on the Duffing oscillator equation shows that the 
nonlinearity, and hence the bifurcation phenomenon, is solely dispersive. 

PACS numbers: 74.70.Ad, 74.25.Fy, 74.25. Qt, 74.25.Jb 



Amplifying in the microwave frequency domain signals 
whose energy are at or close to the quantum limit con- 
stitute an experimental challenge whose pursuit is justi- 
fied both by the fundamental understanding of amplifi- 
cation mechanisms in general [l| and by applications in 
radioastronomy 2] and quantum information processing 
Q. One practical, minimal noise amplification mech- 
anism involves pumping a purely dispersive non-linear 
medium in the vicinity of a bifurcation between two dy- 
namical states. The weak input signals can then be de- 
tected either in the non-hysteretic regime for continuous, 
phase preserving operation or in the hysteresis regime 
where the medium provides latching by switching from 
one state to the other. 

In this Letter, we demonstrate that a microwave su- 
perconducting cavity incorporating a Josephson junction 
can display a bistable regime of operation suitable for 
ultra-low noise amplification. The bifurcation associated 
with this bistability involves a dispersive non-linear dy- 
namical evolution of the fields of the cavity on time scales 
given by its quality factor rather than a slow change in 
its parameters due to dissipation induced heating of the 
cavity material. Importantly, we show that the strength 
of the non-linearity depends on the combination of sim- 
ple electrical characteristics of the junction and the cav- 
ity, both being entirely controllable by fabrication. The 
abscence of dissipation inside the cavity ensures that a 
minimal number of modes are involved in the amplifica- 
tion process, a necessary condition for amplification at 
the quantum limit [1] 

A Josephson tunnel junction constitutes the only 
radio-frequency (RF) electrical element which is both 
non-dissipative and non-linear at low temperatures and 
whose characteristics are engineerable in a wide parame- 
ter range. This property makes it unique for applications 
in ultra-low noise amplifiers, particle detectors, and more 
recently, quantum information processing. The pioneer- 
ing work on squeezing of Yurke et. al. in the 80's [H 



had already showed that the Josephson junction could 
achieve for microwave radiation in the quantum regime 
what non-linear media could achieve for quantum optics. 
Since then, recent works in which a Josephson circuit 
is coupled to a superconducting microwave cavity have 
demonstrated that the strong-coupling regime of cavity 
QED in atomic physics could be attained [6(. Here, we 
are showing that the cavity can become so non-linear that 
the analog of dispersive optical bistability with atomic 
ensembles 0] can be observed. 

In Fig. [1^, we show our superconducting Fabry-Perot 
cavity constructed from a coplanar waveguide transmis- 
sion line by interrupting the conductor with two in-plane 
capacitors playing the role of partially transparent mir- 
rors (see Fig.[Tp, d). The geometric distance between the 
capacitors determines the lowest resonance frequency of 
this resonator. In our experiments, by varying the length 
of the transmission line, the frequency is tuned between 
1.5 and 10 GHz. The quality factor Q of the cavity was 
in the 400-2500 range and was entirely determined by the 
loading through the coupling capacitances. In a separate 
experiment [8[ we have demonstrated that the contribu- 
tion of intrinsic dissipation to Q for this type of resonator 
at our working temperatures (T ~ 300 mK) were negligi- 
ble. In the middle of the cavity, the central conductor of 
the transmission line is interrupted with a gap into which 
a Al-AlOx-Al Josephson junction is inserted (see Fig.[Tfc), 
ensuring maximum coupling between the junction and 
odd cavity modes. Fabrication of the junction involves 
standard e-beam lithography, shadow mask evaporation 
and lift-off. Good metallic contacts between the leads of 
the junction and the Nb resonator were obtained with ion 
beam cleaning prior to the e-beam evaporation of alu- 
minium films. In the following, we will concentrate on 
the lowest order mode for which the analysis developed 
in [11] applies. 

The circuit diagram of our experiment is given in 
Fig. [TJj where the two sections of loss-less coaxial 
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FIG. 1: (a) "Dark field" photographs of a chip sample showing 
the coplanar waveguide Nb resonator implementing a super- 
conducting Fabry-Perot cavity on chip. Asymmetric finger 
capacitors (b)-(d) play the role of partially reflective mir- 
rors. Non-dissipative non-linearity is provided a shadow- 
mask evaporated Josephson junction (c). (e) Schematic of 
the microwave circuit of the experiment. The dashed box 
surrounds the parts which are on-chip. Viewed from the 
Josephson junction the RF biasing circuit appears (f) as an 
impedance Z(u) in series with a voltage source, (g) In the 
vicinity of the first mode of the resonator, the electromag- 
netic environment of the junction can be modeled by an LCR 
combination. 
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FIG. 2: Measured output power (top) and phase (bottom) of 
transmitted signal as a function of increasing drive frequency 
for a range of incoming powers (colors) and for samples 2 (left) 
and 3 (right). The black line corresponds to the critical power, 
separating the single state solution regime from the bistable 
one. The transition between states appears as a vertical jump. 
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FIG. 3: Theoretical predictions corresponding to the protocol 
of Fig. 2 and calculated for the Duffing model, plotted as a 
function of reduced drive frequency 2Q(uid — u>o)/ojo- 



transmission lines represent the coplanar waveguide on 
each side of the junction. The transmission lines are 
parametrized by a characteristic impedance [Zq ~ 50 f2) 
and a frequency loq corresponding to their A/2 (funda- 
mental) resonance. The small input capacitor Ci n cou- 
ples in a source signal from a microwave generator while 
the large output capacitor C ou t couples out the transmit- 
ted signal to the first stage of our amplifier chain. A fac- 
tor of 20 to 35 in capacitor asymmetry ensures that all the 
power escaping from the resonator is collected on the am- 
plifier side. This maximizes the signal to noise ratio of the 
experiment since it prevents the input port to participate 
in the dissipation. Applying the Norton-Thevenin theo- 
rem, our microwave circuit can be seen from the junction 
as an ideal RF voltage generator connected to a series 
combination of the junction and a finite linear embed- 
ding impedance Z (u>) (see Fig. [Tf). For drive frequen- 
cies in the vicinity of the fundamental mode the system 
is very well described by the model shown on Fig. [Tg, 
where Z (uj) is now a series LCR circuit. When com- 
bined with the junction non-linear inductance governed 
by the parameter Lj = (fio/Io, where (fo = h/2e, this cir- 
cuit yields an equation, analogous to the Bloch-Maxwell 
equations described in ref 9] , and which can be analyzed 
in detail [ill ] . In the asymmetric mirrors limit parameters 
of the circuit from Fig. [I] are given by V = Z Lu Ci n Vd, 
R = ZZRlojZC^, L = -. z ° , C = iy 1 m where 
Rl is the 50 f2 load resistance. Note that loq absorbs a 
slight renormalization by the coupling capacitors. The 
quality factor is given by Q — 07 J* . We will come 
back to this model after examining the data. 

The basic principle of our experiment is to measure 
the transmission characteristics of the cavity as a func- 
tion of both frequency and power. We studied three dif- 
ferent chip samples for which the relevant parameters 
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are given in Table HI They were anchored to the final 
T=300 mK stage of a 3 He refrigerator. While we have 
mainly conducted CW measurements, we have also con- 
ducted measurements in which, keeping the frequency 
fixed, the power is ramped sufficiently rapidly to probe 
the internal dynamics of the cavity. Since we would like 
to demonstrate that the cavity is strongly non-linear even 
for small radiative energies stored inside it, the challenge 
of the experiment is to perform precise transmission mea- 
surements at very low power. A vector network analyzer 
sends a CW microwave signal on the small capacitor side 
and analyzes the transmitted signal coming out from the 
large capacitor side, after they have passed through the 
amplifier chain which includes one cryogenic HEMT am- 
plifier at 4 K and two circulators placed at 300 mK and 
4 K. 

In Fig. [5] we show the typical transmitted signal ampli- 
tude (top panels) and phase (bottom panels) as a func- 
tion of drive frequency for different input drive power 
and for two resonator frequencies (samples 2 and 3). 
The drive frequency was swept in 800 ms. At low drive 
power Pi n — 2TTZ^d the response of the resonator is a 
Lorentzian centered on the resonator frequency. As drive 
power is increased, the resonant response sharpens on 
the low frequency side and softens on the high frequency 
side. At a critical drive power Pi n = P c , the maximum 
slope on the low frequency side diverges. Upon increas- 
ing the drive power further, a abrupt transition develops, 
whose position shifts to lower frequency with increas- 
ing power. The two panels of Fig. [2] demonstrate that 
the phenomenon presents itself in identical manners for 
two resonators despite their different resonant frequencies 
and quality factors. For clarity the response at the crit- 
ical power is highlighted in black. At even higher drive 
power (data not sho wn) , the system becomes chaotic as 
described elsewhere [id ]. 

These results should be fully governed by the static 



solutions of the Duffing oscillator as described in 
Indeed, simple analytical relations for the transmitted 
power and phase at different drive powers show excellent 
qualitative agreement with the measurements as shown 
in Fig. [3] We now focus on the locations in the (u^, Pi n ) 
plane where the susceptibility is maximum, i.e. loca- 
tions where dP ou t/duJd is maximum. When P m > P c , 
these should correspond with the bifurcation points. We 
will also follow the point of maximum P ou t- These re- 
sults are plotted for the three samples and three dif- 
ferent values of the DC critical current of the Joseph- 
son junctions patterned in our samples. We plot them 
against the reduced frequency f2 = 2Q(ujq — uJd)/uio and 
Pin/Pc, both obtained independently for each sample. 
As expected, the various experimental points fall on a 
universal theoretical curve when properly scaled. Based 
on (llT | we expect the bifurcation power to be given by 







sample 1 


sample 2 


sample 3 


measured 


wo (GHz) 

Q 

C in (pF) 
Pc(dBm) 


1.828 

2360 
2.3 ±0.2 
-85.5 ± 1 


1.8177 
1500 
1.5 ±0.5 
-76.5 ± 1 


9.246 
460 
1.0 ±0.1 
-70 ± 1.5 


inferred 


L (nH) 
C (pF) 
R (fi) 

h (M) 


6.7 
1100 
0.034 
1.6 ±0.2 


6.8 
1100 
0.053 
1.5 ±0.5 


1.3 
200 
0.17 
2.9 ±0.35 


expected 


Jo (juA) 


1.3 ±0.1 


1.6 ±0.1 


3.35 ±0.15 



TABLE I: Table showing the measured values of parameters 
for samples studied in experiment, as well as inferred and ex- 
pected parameters. d n is obtained either from the measure- 
ment of other identical resonators with no Josephson junc- 
tions (samples 1 and 2) or by measuring the ratio of input 
and output powers leading to criticality and assuming that 
the quality factor is set by C ou t (sample 3). 
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red lines respectively) and the power of highest derivative 
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FIG. 4: Logarithmic plot of the highest derivative dP ou t/du>d 
as a function of reduced parameters Q and Pi n jP c . Lines de- 
pict theoretical prediction for the Duffing model. Triangles, 
circles and crosses represent respectively the measured bifur- 
cation power, its highest derivative (below P c ), and the data 
of the hysteresis measurement. Inset: derivative of the out- 
put power with respect to drive frequency, normalized to its 
maximum as a function of drive frequency and power. Note 
the white line delimits the point where the output power is 
maximum (the derivative changes sign). 



below P c to be Phd/Pc = - \ (dashed line). They 
are also shown on a plot of the output power in the inset 
of Fig. [4] There, dP ou t/dujd is normalized to its maxi- 
mum and plotted as a function of the drive frequency and 
power. The maximum output power (below P c ) is shown 
as a white line defined by P m ax/Pc — ^75^ an d coincides 
with the change of sign of dP ou t/duid- Note that these 
measurements do not probe the hysteresis since the fre- 
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quency is swept only in the forward direction. To verify 
the hysteretic behavior of the phenomenon, we used tri- 
angular power sweeps for several frequencies. We were 
able to probe the power and frequency dependence of 
both bifurcation points. The resulting data is shown as 
stars on Fig. 2J The deviation from predictions is most 
likely due to the proximity of the lower bifurcation cur- 
rent amplitude to Iq, situation not well addressed by the 
Duffing model. Note that the upper and lower bifurcation 
current in this RF experiment correspond to the switch- 
ing and retrapping current, respectively, in DC Josephson 
IV measurement. 

We now discuss how, from the measured value P c 
of the critical power, we can infer Iq. In contrast to 
a DC experiment, where one can characterize the bi- 
asing circuit with great precision, an absolute calibra- 
tion of the environmental impedance Z(lu) at RF fre- 
quencies is arduous. From the values of Ci„, Zq, L 
and Q obtained from the low power (linear) measure- 
ments, as well as the relations crit — ancl 
V ° = ^( i E3Q L ) 3/: Wo, we extract 7 
the measured value of P, 



V dcrit/ 2R L- 



ZoCinUJQ 

ipo/Lj from 
In table 1 we 



show that Iq obtained in this way is consistent with the 
normal resistance test of junctions fabricated in the same 
batch. We can thus verify that the Duffing model hypoth- 
esis I c /I = -|=y^±^ii < 1 for the RF critical current 

amplitude I c is well satisfied (llj . 

Other experiments on an AC biased Josephson junc- 
tion have been realized 

HEBE!, albeit with 

a lower 

degree of control over the electrodynamic environment. 
Nonlinear effects involving superconduct ing weak links 
in resonators have been reported 15, 3, 17 1. Moreover, 
the Duffing oscillator physics has been shown to appear 
in several other systems including nanomechanical res- 
onators [Hj and relativistic electrons However, our 
realization offers an unprecedented opportunity to study 
dynamical driven systems in the regime where quantum 
fluctuations are dominant 0, [2(|. This regime is yet 
unexplored experimentally. 

The excellent overall agreement between experimental 
results and theoretical predictions allows us to eliminate 
non-linear dissipative effects as the cause for the bifur- 
cation. For instance, if the dissipation in the resonator 
would increase with power the resonance would become 
broader, its maximum would decrease and the resonance 
curves will eventually cross each other, an effect never 
observed in our experiment. 

The main application we envisage in the short term 
for this phenomenon is parametric amplification at mi- 
crowave frequency which would reach the quantum limit 
at the practical level, not simply at the proof-of-principlc 
level. It could also used as a superconducting quantum 
2l|. Because the cavity frequency is mainly 



domain several cavities and multiplex the measurement 
or readout of a collection of entities. 

In conclusion, we have observed the dynamical bifur- 
cation of a superconducting microwave cavity incorpo- 
rating a non-linear element in the form a Josephson tun- 
nel junction. This bifurcation is analogous to the optical 
bistability observed with atoms in a Fabry-Perot in QED 
experiments. Our system is therefore at the crossroads 
between the physics of dynamical systems and quantum 
optics. Our level of control on the various parameters 
of this bistable phenomenon opens the door to amplifi- 
cation at the quantum limit, squeezing and other analog 
quantum information processing functions in the strong 
and ultra-strong coupling regimes. 

The authors are grateful to Daniel Esteve, Denis Vion 
and Steve Girvin for helpful discussions. This work was 
supported by NSA through ARO Grant No. W911NF- 
05-1-0365, the Keck Foundation, and the NSF through 
Grant No. DMR-0325580. 
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